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1 . 


IDENTIFICATION 


1.1 Number: 0CDA-8-08 (P-08-01) 

1.2 Title: DISTANCE AMD BEARING 

1.3 Date: May 1968 

1.4 Computer: PDP-8 

1.5 Language: PAL III 

2. ABSTRACT 

The distance and bearing are calculated for two positions 
specified by their latitudes and longitudes. 

3. REQUIREMENTS 

3.1 storage: 5 - 7, 40 - 104, 200 - 577, 4600 - 7577 

3.2 Subprograms and/or Subroutines 

3-2.1 F.P. Package (Digital 8-5D-S) 

3-3 Equipment: Teletype 

4. USAGE 

4.1 Loading: The Binary loader (Digital-2-U) is used to load the F.P. 
Package followed by the program. 

4.2 Calling Sequence: N/A 

4.3 Switch Settings: N/A 

4.4 Start up and/or Entry: The program is started at loc. 200^. Then the 

positions are typed on the teletype. . 

4.5 Errors in Usage and Recovery: A space is the only legal terminating 
character. If any other terminating character is read, the preceding 
number is ignored, (i.e. all digits typed between the last space and 
the illegal terminating character). To continue it is necessary to 
retype the number. 

5. RESTRICTIONS N/A 

6. DESCRIPTION 

6.1 Discussion: The program calculates the distance and bearing between 
positions. Initially the program requires two positions, but only 
one thereafter. The distance and bearing between the last two 
stations typed are calculated. All I/O is on the teletype. 

6.2 Examples and/or Applications: Given in Fig. 1 

7. METHODS 


7.1 


Discussion: N/A 





44 23.8 63 28.2 
44 15.5 63 19.5 
43 52.8 62 53.2 
43 28.8 62 27*2 
43 10.6 62 06.3 
42 50.8 61 44.1 


+ 0019.24 + 036.96 
+ 0054.78 + 039.87 
+ 0056.53 + 038.18 
+ 0043.98 + 039.98 
+ 0047*48 + 039.45 


Figure 1 
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Algorithm: The following equations are used to find the geodetic 
distance between two points given their latitudes and longitudes: 

= aA ci> m _ ,, 

A. 

S co-^oC = - = \/ 

- (I - efW cp m f 
* 

R , ( i- 

^ r*q- 

°-(\ - e 1 -) a^ry. X* 


where X, ^ X z : 

tP.A’ 

S: 

<A,= 

oC • 


e : 
a : 
b : 

^ I" : 


longitude of the two points 
latitudes of the two points 
distance bwtween positions 
mean latitude 
azimuth 

2 

eccentricity e = .0067686275 
equatorial semi-axis 
polar semi-axis 

permits to be expressed in seconds. 


A A = \-\ l _ 

& $ = 4> - 


i _= .032559 

a sin I" 


1 _ = .032559 

a (l-e^) sin I" 

8. FORMAT 

8.1 Input Data: Each position is given by four numbers which represent 
the latitude and longitude in degrees and minutes. The terminating 
character for each number is a space. 

8.2 Core Data: N/A 

8.3 Output Data: The O/T record consists of the distance in kilometers 
and the bearing in degrees (+xxxx.xx +xxx.xxJ ). Note that if the 
bearing is a multiple of 90, a series of x's will be printed. 

9. EXECUTION TIME N/A 
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PROGRAM 


10.1 Core Map: 



10.2 Dimension List: N/A 

10.3 Macro, Parameter, and Variable List: This is included with the 
listing of the program. 

10.4 Program Listing: Attached at the end of the program write-up. 

11. DIAGRAMS 

11.1 Flow Chart: Shown in Fig. 2 

12. REFERENCES 

12.1 Other Library Program: 

12.1.1 F.P. Package (Digital 8-5D-S) 

12.2 Digital Manuals: 

12.2.1 F.P. Manual (Digital 8-5-S) 

12.3 DECUS Programs: 

12.3.1 Modification to Fixed Point O/T in the PDP-8 F.P. Package 
(Decus - No. 8-44) 

12.4 Textbooks: 

12.4.1 U.S. Coast and Geodetic Special Publication #8 
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